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Recent progress in nanotechnology allows to engineer hybrid mesoscopic devices comprising on 
chip an artificial atom or quantum dot, capacitively coupled to a microwave (superconducting) 
resonator. These systems can then contribute to explore non-equilibrium quantum impurity physics 
with light and matter, by increasing the input power on the cavity and the bias voltage across the 
mesoscopic system. Here, we build such a prototype system where the artificial atom is a graphene 
double quantum dot (DQD). Controlling the coupling of the photon field and the charge states 
of the DQD, we measure the microwave reflection spectrum of the resonator. When the DQD is 
at the charge degeneracy points, experimental results are consistent with a Kondo impurity model 
entangling charge, spin and orbital degrees of freedom. The light reveals the formation of the Kondo 
or Abrikosov-Suhl resonance at low temperatures. We then study the complete phase diagram as a 
function of gate voltages, bias voltage and microwave input power. 


I. INTRODUCTION 

Circuit quantum electrodynamics (cQED) with quan¬ 
tum dot geometries has gained rapid development re¬ 
cently [IHil offering a platform to address interaction 
effects between photons and electrons on-chip. These 
settings typically involve an electronic nano-circuit, such 
as a single or double quantum dot (DQD) coupled to 
source/drain leads and to an electromagnetic resonator. 
A bias applied across the nanocircuit can cause a current 
flowing as well as dissipation, which will result in chal¬ 
lenging physics when coupled to a photon field. These 
hybrid systems are not only of great practical relevance 
for quantum information technology | 2 l[T 2 l [23 but also of 
fundamental interest as platforms to explore exotic quan¬ 
tum impurity physics with light and matter [26) . In this 
respect DQD-resonator devices are particularly interest¬ 
ing, as they allow deep quantum limit transport inves¬ 
tigations, including electronic entanglement [smiiiiHi- 
[2D], DQD maser nni m na, photon emission statistics 
Ha heat engines m- From a more fundamen¬ 
tal perspective, a DQD may host interesting degener¬ 
ate states, entangling charge, spin and orbital degrees 
of freedom, which can turn through quantum fluctua¬ 
tions into exotic many body states [SZHSa- We note re¬ 
cent progress concerning the strong coupling of spin to 
photons through artificial spin-orbit interactions [la or 
polarization rotation m- Out of equilibrium charge dy¬ 
namics in a hybrid cQED architecture na, and photon 
emission from a cavity-coupled DQD [laiia have been 
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reported. In this paper we present evidence for exotic 
Kondo correlations in the microwave response of a DQD- 
cavity coupled to leads, a unique signature of many-body 
quantum impurity physics with light, which has not yet 
been achieved in these hybrid devices. 

We probe the response of a DQD-resonator device 
by measuring phase and amplitude of the reflected mi¬ 
crowave signal as a function of bias and gate voltages 
as well as the response of the resonator as a function of 
the input microwave power. Our data reveal how unique 
features of the electronic transport appear in the pho¬ 
ton spectroscopy. In particular at the charge degeneracy 
points of the DQD, the phase response of the reflected 
light signal in this geometry shows a robust tt phase shift 
upon increasing the bias voltage up to a threshold set by 
the charging energy, while the amplitude response shows 
additional satellite peaks at low bias. We show that the 
experimental findings are in qualitative agreement with 
an effective light-matter quantum impurity model in rela¬ 
tion with Kondo physics [IS]. More precisely, we extend 
the Kondo model in Ref. m applied to the DQD at the 
charge degeneracy points to the case of graphene DQD. 
SU(4) Kondo effect was recently observed in a GaAs 
double quantum dot |33|. We justify that the decou¬ 
pling of the antisymmetric wave-function in terms of the 
two graphene sub-lattices leads to a fully screened SU(4) 
Kondo model in the low-temperature limit. We will show 
how the experimental datas are consistent with this inter¬ 
pretation at very low temperatures. We shall also discuss 
temperature effects. We will also assume that the chem¬ 
ical potential does not lie at the charge neutrality point 
in the graphene leads, and justify this assumption. For 
an introduction on Kondo physics, see Refs. |45H5T|; for 
Kondo physics in quantum dots, consult Refs. 
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FIG. 1. (color online), (a) Micrograph of the DQD gate 
structure, (b) Sample structure of a typical etched graphene 
DQD. The dc voltages used to control the charge numbers 
in the DQD are applied via left and right plunger (LP and 
RP) gates. A quantum point contact with a source (Sq) 
and drain (Dq) channel and a tuning gate (Q) is integrated 
near the DQD. (c) Circuit schematic of the hybrid device. 
The half-wavelength reflection line resonator is connected to 
DQD’s left dot (LD) at one end of its two stripelines. The 
right dot (RD) is connected to the drain. A microwave signal 
is applied to the other end of the resonator, and the reflected 
signal is detected using a network analyzer. 


To our knowledge, our paper constitutes the first report 
of the realization of Kondo physics with graphene dou¬ 
ble dots at the charge degeneracy point exhibiting charge 
and spin entanglement m- Previous works on graphene 
double quantum dots were carried both in series shape 
m and in parallel [SB]. The Kondo effect has also been 
observed in graphene with point defects (vacancies) [55] . 
Similar to a three-lead or four-lead geometry O EU] , the 
cavity which is weakly-coupled to the mesoscopic system, 
allows to reveal certain aspects of the spectral function 
on the DQD. Josephson-Kondo circuits have also been 
proposed to probe Kondo physics of light [BlHHS] . We 
note related progress from the experimental point of view 
both in cQED [511 - 155] and in mesoscopic systems m- 

Our device (shown in Fig. 1) is mounted in a dry di¬ 
lution refrigerator, its base temperature is about 30 mK. 
Two DQDs, made of few-layer etched graphene, are cou¬ 
pled to the resonator through their sources El, how¬ 
ever, only one of them is used in this experiment while 
the other is grounded all over the experiment. A vector 
network analyzer (VNA) is used to apply coherent mi¬ 
crowave driving tone and measure the reflection signal 
(see Fig. 1(c)). 

The DQD may have several effects on the resonator: 
dispersive effects 0191 can cause a shift on the resonance 
frequency as well as dissipation effects [5]. Derived from 
the input-output theory [miia, the reflection coefficient 


FIG. 2. (color online) (a) Gharge stability diagram of the 
DQD measured by the dispersive readout using the resonator 
and (b) Goulomb diamond measured by transport experiment 
and revealing the charging energy. The gate lever arms are 
~ 10 %. 
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FIG. 3. (color online) Results for point A. (a) Best fit of the 
phase as a function of the frequency when DQD bias voltage 
I4d=l mV, from which the resonance frequency ujq, Ki and Ke 
can be obtained, (b) Amplitude response as a function of the 
driving frequency, for various bias voltages (0, 0.5, 0.6, 0.7, 
0.8, 1 mV). The red dashed line shows the fitted resonance 
frequency, shifted by the electron transport, (c) Ki and Ke as 
functions of the DQD bias voltage, (d) Total dissipation and 
load quality factor Ql = ojQ/{Ki + Ke) as a function of the 
bias voltage. 


of the resonator is defined as [21 [73]: 


S'ii(a;) 




( 1 ) 


where oj is the driving frequency, oJq is the cavity res¬ 
onance frequency (taking into account the shift from 
the coupling with the mesoscopic system), Ke represents 
the dissipation mechanisms from the light-matter cou¬ 
pling and Ki incorporates the dissipation from the pho¬ 
ton system only. In the present experiment, these pa¬ 
rameters can be explicitly tuned with the bias voltage 
applied across the mesoscopic system and a precise anal¬ 
ysis will be required to fit the experimental datas. The 
total dissipation rate on the cavity is k = Ki + Ke- 
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FIG. 4. (color online) Results for point A, for the ampli¬ 
tude of the reflection coefficient 5ii(a;o) as a function of the 
source-drain bias; the zoom focuses on low-energy Kondo fea¬ 
tures. Theory results based on an effective quantum impurity 
model. We take Tk = 550mK (extracted from the temper¬ 
ature analysis) and Ec = OAmeV from Fig. 2, uq = GGHz, 
A/tuo = 2.10“^, a = 3.10“^ and to = Tk- For |Vsd| > Tk-, 
we model decoherence effects on the Kondo resonance with 
a (decoherence) rate r(I4d) which grows rapidly with 114^1 
|68H70| and saturates around Ec- 


diamond in Fig. 2(b) is obtained from transport measure¬ 
ments. 

In Figure 3, we plot the amplitude and phase response 
as a function of the driving frequency at point A. Using 
a fit with Eq. (1), we can deduce Wg, Kq and Ke [HI [Z3] 
and their evolution with respect to the bias voltage. Fit¬ 
ting the Usd=0 curve, we obtain a;Q(I4d = 0)=6.35080 
GHz, Ki(V;d = 0)=0.498 MHz and = 0)=1.353 

MHz. For t4d=l mV curve, we obtain a;o(I4d = 
lmU)=6.35084 GHz, Ki(V;d = 1 toU)=1.278 MHz and 
Ke(V 3 d = 1 toU)=0.962 MHz (Fig. 3(a)). Gompared to 
the Usd=0 curve, there is a 40 KHz frequency shift, which 
is caused by the nonlinearities induced by the nonlinear 
electronic transport [ 2 ( 1 ] . Figure 3(c) shows the compar¬ 
ison of Ki and Ke for various bias voltages. From zero 
bias to Vsd=lmV, k; increases by 0.78 MHz and Kg de¬ 
creases by 0.391 MHz. In this description, k; and Ke do 
not depend on the frequency. 


II. THE CHARGE DEGENERACY POINT 
A. Experimental Results 



FIG. 5. (color online) Experimental and Theory results for 
the phase of S'ii(cdo) at Point A. The phase behavior (roughly 
TT or 0) can be understood from the evolution of Ke and Kt with 
the bias voltage. 


Sii can be measured by its amplitude (A = |S'ii|) and 
phase {(j) = arg{Sii)) components or by the imaginary 
(5 to(S'ii)) and real (3fie(S'ii)) parts through the VNA. 
We then set the DQD to different states, to see the ef¬ 
fect of the DQD on the resonator. We can first obtain 
a charge stability diagram of the DQD by a dispersive 
readout measurement [51ini[see Fig. 2(a)]. The Coulomb 


Now we discuss additional experimental results for 
point A, at the resonance frequency (w = Wg). Figures 
4 and 5 show the evolution of the amplitude and phase 
response at the resonance frequency with respect to the 
bias voltage. 

A phase shift Acj) ~ tt is observed from bias voltage 
I4d 0 to Usd ~ ±0.4 mV. The behavior of the phase 
response can be understood from very general arguments 
and Eq. 0 - Indeed at small bias voltages, one expects 
that the coupling of the cavity resonator with the 
electronic system is large then producing Ke ^ Ki and 
a phase of tt in the reflected microwave signal. On the 
other hand, at large bias voltages, the current through 
the DQD should be large such that the dynamics of 
the electronic degrees of freedom is much faster than 
the dynamics of the photon resulting in Ke ^ Ki and a 
phase which is zero (or 27r) |^. This relatively simple 
argument allows to explain the phase shift of — tt (tt) 
in the reflected microwave signal when increasing the 
bias voltage across the mesoscopic electron system, and 
the datas of Figure 3(c) confirm this evolution. The 
amplitude (in dB) shows also two pronounced dips 
around Usd ~ ±0.4 mV, which is consistent with this 
argument. 

It is important to note two points concerning these 
measurements. First, the absolute value of the voltage 
corresponding to the phase shift and the amplitude dips 
corresponds to the the charging energy Ee ~ 0.4meU of 
the DQD measured through dc transport and leading to 
the Goulomb diamond (see Fig. 2(b)). Interestingly, we 
also notice the low bias regime shows a (slightly asym¬ 
metric) dip structure (see Fig. 4(b)). As explained be¬ 
low, we attribute the low bias dip structure as well as 
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FIG. 6. (color online) Upper pannel: temperature evolution 
of the differential conductance at point A, showing two peaks 
at finite voltage Vsd ~ ±0.05mU and a zero bias dip cor¬ 
responding to a small residual conductance, which smoothly 
disappears above a crossover Kondo temperature of about 
Tk ~ 550mK. Curves are displaced from the base ref¬ 
erence by a constant offset. For an SU{4:) Kondo model, 
the zero-bias anomaly in the conductance does not occur at 
Vsd —>■ 0, but rather at voltages in relation with Kondo physics 
|29| . The Kondo anomalies disappear at a temperature scale 
Tk ~ 550mK. In addition, when Vsd —^ 0, at low temper¬ 
atures, one predicts Go = dI/dV{Vsd = 0) ~ {t'/TK)'^2e^/h 
|27| which confirms that t' ~ 0.1 — O.STk <C Tk- Lower pan- 
nels: Current (a) and differential conductance (b) at point B 
at r=30mK, showing no zero-bias anomaly. 



FIG. 7. (color online) Temperature dependence of the mi¬ 
crowave reflected signal amplitude at point A. Curves are 
displaced from the base reference by a constant offset. (We 
emphasize that in this cooling procedure, the effective elec¬ 
trical temperature of the sample is around lOOmA whereas 
the dilution refrigerator temperature T can reach lower tem¬ 
peratures. The datas do not change much on the figure for 
temperature T smaller than lOOmA.) 



the robustness of the tt phase (at low temperatures) to 
Kondo physics. A phase of tt = 2 x 7r/2 is expected in 
the reflected light wave, associated with a bound state 
formation [43]. This phase of tt is not in disagreement 
with the phase of 7r/2 in the transmission predicted in 
Ref. [25]. 

For completeness, we also present in Fig. 6 the trans¬ 
port datas at the points A and B (and the evolution 
of the transport datas with respect to temperature for 
the case of point A). At the charge degeneracy point A, 
these measurements show at low temperatures a very 
rich small-bias structure, well inside the charging energy 
bands, with two peaks at finite voltage Vsd ~ iO.OSmU 
and a zero bias dip corresponding to a finite small 
residual conductance. Upon heating these features 
smoothly disappear above a crossover temperature of 
about 550mK. In Fig. 7 and 8, we show the tem¬ 
perature evolution of the amplitude and phase of the 
reflected signal. In this cooling procedure, the effective 
electrical sample temperature is around 100m AT whereas 
the sample environment temperature of the dilution 


FIG. 8. (color online) Evolution of the phase of the microwave 
reflected signal at point A. The cooling procedure is similar 
to that used in Fig. 7 for the amplitude. Interestingly, one 
observes that the phase converges to tt for T < Tk = 550mK 
for Vsd —>■ 0 and that the phase evolution is sensitive to the 
Kondo temperature. In addition, experimentally, for inter¬ 
mediate temperatures (we remind that the effective electrical 
temperature of the sample is around 100mA in this cooling 
procedure), the phase smoothly drops to zero for |I4tj| > Tk- 


refrigerator can reach lower temperatures. Therefore, 
the datas of Fig. 6, Fig. 7 and Fig. 8 do not evolve much 
for temperatures of the refrigerator below lOOmAT. 


B. Exotic Spin and Orbital Kondo model 

To support the experimental results, we provide a the¬ 
oretical model to describe the precise evolutions of Ke 
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and Ki with the bias voltage. 

At the point A, the two low lying charging states (1, 0) 
and (0,1) — with one electron either on the left or on 
the right dot — are degenerate since it costs the same 
energy to have one extra electron on the left or on the 
right dot [74] . Other configurations on the DQD cost 
higher energy. The low-energy theory can then be built 
by introducing, in addition to the spin = ±1/2 of 
the electron delocalized on the DQD, the orbital pseu¬ 
dospin quantum number = ±1, projecting on these 
two allowed states which are further coupled by a finite 
tunneling term t'T^. These four quantum states of the 
DQD generate exotic Kondo physics when coupled to the 
two channels conduction electrons. 

As a result of second-order tunneling processes, the 
electron leads are coupled through the Hamiltonian Hi = 
Hkin ± Hkondo ± ^assist ± ^orbital [^'7] ■ Here, Hkin rep¬ 
resents the kinetic energy in the two leads. Second-order 
tunneling processes are classified through purely Kondo 
terms involving spin flips, orbital contributions changing 
the lead index from say I = L (left) to I = R (right) 
and flipping the charge state on the DQD, and assisted 
tunneling processes entangling the charge and spin HZj: 

HKondo ± Hassist = ( 2 ) 

+ Q± {t^T~ aij) ± h.c^ 

HorUtai = ^ + y± + h.c.)) , 

where tpoi is the Fourier transform of Cka-i which refers 
to an electron operator with spin a in lead 1. Here, S 
corresponds to the spin operator of the electron delocal¬ 
ized on the DQD, a and r act on the spin space and 
orbital (lead) space of the electron operators. The cou¬ 
plings above can be written in terms of the tunneling rate 
to the right/left lead r± and the charging energy of the 
DQD Ec respectively as [27] J = Qz = + r_)/4ifc, 

Qi = 14 = yD/rl/Fic. 

The model presented here from Ref. [27| can be ex¬ 
plicitly refined to incorporate the precise geometry of 
the graphene lattice (leads) [75] . Let us consider the two 
sub-lattices of the honeycomb lattice A and B. The 
above analysis remains valid assuming that the symmet¬ 
ric wave-function tpoi = ± '4’Bai) couples to 

the double dot through second order tunneling processes, 
whereas the anti-symmetric component would be discon¬ 
nected from the double dot and therefore can be ignored. 

The kinetic term leading to a Dirac type linear dis¬ 
persion can also be diagonalized in this basis if we as¬ 
sume a symmetric configuration for the A and B sites 
with respect to the direction of the tunneling. The use 
of either open or periodic boundary conditions in the or¬ 
thogonal direction then lead to a one-dimensional kinetic 
Hamiltonian along the direction of the tunneling, which 
is diagonal in this basis [76l [77]. We argue that a small 
geometrical asymmetry is unimportant, as the symmetric 
term would dominate in the renormalization procedure. 


Below, we assume that the density of states in graphene 
is finite, meaning that the chemical potential does not lie 
at the charge neutrality point where the density of states 
vanishes (The Kondo temperature would vanish, which 
does not seem in agreement with experimental results). 

As shown in Ref. m and in related models [28ll30| . 
this DQD model at point A flows to a Kondo-type fixed 
point with a Fermi liquid ground state [75H5D] . which is 
quite robust to (charge) noise effects [ST]. Experiments 
have reported the possible occurrence of an emergent 
SU{4) symmetry at the low-energy fixed point, in simi¬ 
lar geometries, due to the entanglement between spin and 
orbital degrees of freedom [STTfSS] . The remarkable fact 
is that the bare parameters r_|_ and F_ are replaced by 
a single parameter describing the low-energy fixed point, 
namely the Kondo energy scale Tk- In this description 
the spectral function of the DQD can be modeled by an 
effective resonant level model [29], associated with spin¬ 
less fermionic operators d and db These operators d and 
rfl are defined in correspondance with the processes in¬ 
volving the orbital operators T~ and respectively. 
This spinless description is legitimate as the two spin 
components add up in the conductance, and the cavity 
is only sensitive to the orbital degree of freedom. In this 
strong coupling regime, the system can then be treated 
by analogy with the case of one single dot studied in Ref. 
[7S1 the only difference being that the position eg and 
the width F of the electronic level are determined by the 
SU(4) fixed point. At small bias voltages Vsd < Tk, the 
Kondo resonance is described by the spectral function 

m- 

where eg is of the order of Tk [IHl [50] . We stress that 
even though we describe the effective resonant level with 
spinless operators, the spin of the electrons is a highly 
relevant quantity: the joint effect of the orbital and spin 
degrees of freedom lead to a Kondo resonance at eg of 
the order of Tk, above the Fermi surface. In contrast, a 
purely spin Kondo effect would yield eg = 0, and trans¬ 
port properties would be very distinct. 

We first remark that the DC transport datas exposed 
in Fig. 6 is in accordance with our model. The transport 
datas for the point A indeed yield a zero-bias anomaly at 
Vsd ^Tk (see the upper pannel of Fig. 6), which is con¬ 
sistent with this form of the spectral functions, with eg 
above the Fermi energy of the reservoir electron leads. 
The conductance across the DQD at small bias volt¬ 
ages essentially follows 2e^/hTKp(bd = e\Vsd\){t' /TkY , 
, where h is the Planck constant and Tk is the Kondo 
energy scale. The spectral function gives the density 
of states accessible for an incoming electron, the pref¬ 
actor {t'/T kY corresponds to the transition probabil¬ 
ity between the two dots. This also leads to a finite 
conductance at Vsd —>■ 0, as observed in Fig. 6. In 
the Kondo limit, one predicts a differential conductance 
Gg = dI/dV{Vsd = 0) ~ (tV74)^2e^/h [57], which leads 
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to the prediction t' ^ 0.1 —O.STff. In addition, the 5'C/(4) 
Kondo theory predicts that the conductance evolves (in¬ 
creases) linearly with the bias voltage for 14 d —>■ 0 [H], as 
found in the experiment (Fig. 6). From the temperature 
analysis of the datas, we also infer that the Kondo energy 
scale Tff, at which the zero-bias peaks in the differential 
conductance (Fig. 6) and the low-bias features disappear 
in the light reflected signal (Fig. 7 and Fig. 8) is roughly 
around 550mK. 

By contrast, the differential conductance at point B 
does not show such anomaly at Vgd ~ as can be 

seen in the two lower pannels of Fig. 6. At point B, the 
levels of the two dots are not degenerate : this energy 
splitting is the analog of an orbital magnetic field along 
the z direction, which suppresses the orbital Kondo effect 
and suppresses the conductance at —>■ 0. 


C. Light-Matter Coupling 


The cavity and light-matter coupling are described 
through the Hamiltonian 


H2 


ojoa^a + A 


_ k(7 


(a -I- o'^). 


(4) 


where we fix the Planck constant h = h/(27r) to unity. 
The cavity is capacitively coupled to the left lead only. 
The key point is that one can re-write the light-matter 
coupling as a quantum dot-cavity coupling of the form 
XT^{a^ + a). To show this, one can apply a unitary trans¬ 


formation U = exp 



with 


I _ x{a a ) Hamiltonian H = Hi + H 2 + t'T^. 

' lOJQ 

The part of the transformed Hamiltonian involving the 
cavity takes the form: 


+ T+e"*'^) - XT^{a + a'f)+ woa^a(5) 


The phase factor entering the tunnel coupling between 
the dots appears in fact at all the investigated points 
in the phase diagram; the total Hamiltonian then reads 
H = Hi + Hcav, see Appendix A. 

We can proceed and compute the reflection coefficient 
using the input-output theory [niEg. The input power 
couples to the x oc (a -|- a'^) variable of the cavity. There¬ 
fore, at small input power, the reflection coefficient reads 

m- 

Sii{uj,Vsd) = -l-h2iJ(a;)x^,,(a;, Wd). (6) 


The susceptibility for the photon as well as the photon 
self-energy are defined as [26] 

^ - UJoUR{uj) + iJ{uj)uJo ’ 

where the photon self-energy H'^ = ab¬ 

sorbs the light-matter coupling. Using these identifica¬ 
tions, for frequencies close to Wg, we recover Eq. (1) with 


Ki = J{ujQ{Vsd)) = 27raa;5 and Ke = 3mn^(wo). At 
zero temperature, the function J(uj) describes explicitly 
the dissipation of the cavity mode due to the coupling 
with the long transmission line which transports the in¬ 
put signal (and we can safely neglect the effect of the 
small number of photons produced in the cavity for volt¬ 
ages larger than Tx due to the current flowing across the 
dot dSj). Theoretically, we assumed the specific form 
J(u!Q{Vsd)) = ^irauiQ which is typical of a transmission 
line (ohmic bath) in the quantum limit [82U88] (with an 
ultraviolet frequency cutoff much larger than Tr-); a is 
a dimensionless parameter proportional to the resistance 
of the transmission line. The coupling between the left 
quantum dot and the cavity produces an extra dissipa¬ 
tion source for the cavity. In the text above, the cavity 
frequency is defined auto-coherently by 

w*(Wd)2 =ujI+ K(Wd)]. (8) 

The measurements of the renormalized resonance fre¬ 
quency carried out on our device match well the theoret¬ 
ical expectations. We note in particular that a;g(Usd) < 
ujQ increases with Vsd — reflecting that SReH^ becomes 
less important — and tends to ujq at large bias [26] . 
as can be seen in Fig. 3(b) . This also implies that 
JIcoq) or Ki becomes more substantial when increasing 
the bias voltage, as experimentally confirmed in Fig. 3(c). 
In addition, from the experimental Fig. 3(b) and for¬ 
mula (7), we estimate roughly the light-matter coupling 
A/wg ~ 2.10“^ which justifies the perturbative analysis 
above. 


D. Effective Quantum Impurity Model description 
of Photon Transport 

Computing the photon self-energy for the exotic spin- 
orbital photonic Kondo model introduced in the previ¬ 
ous section is a challenging task. To proceed we use the 
fact that the electronic sector flows to a Kondo-Fermi- 
liquid fixed point [2HHS0]j where the spectral function 
on the DQD can be modeled by an effective resonant 
level model, as exposed previously. We take a bias de¬ 
pendent width F to include voltage induced decoherence 
effects [SHI El] on the pseudo-fermion d. We then show 
thanks to Keldysh perturbation theory [26] that the pho¬ 
ton self-energy H^ = IfteH^ — i^mH^ is related to the 
density-density response function of a purely electronic 
Anderson Impurity Model. We have at second order in 
the light-matter coupling 

n'^(^) = 2 j ^ - H) -I- Gq{uj)Gq{u} + Q)] , 

(9) 

where Gq, Gq and Gg are respectively the Retarded, 
Advanced and Lesser Green’s functions of the electronic 
level associated to the Kondo resonance. Using this de¬ 
scription of the spectral function [SS] on the DQD we 
evaluate the photon self energy 3mn^(a;Q, Usd) at the 









7 


renormalized resonance frequency Wg (14d) 


3mn^(a;5) = A^/r(a;o) ^ aarctan 

Q:,a=± 


Ma - Co + awg 


+AVrK) XI 


- -u <^0 

,a—± 


(Mo ~ eo + + r 

(Ma - eo)^ + 

( 10 ) 


where /r(w5) = r/(47r2r2 + tt^Wq^), /^a = aT4d/2 and 
r is bias-dependent: F = Tk the Kondo energy scale 
(of the order of Wg) at small biases and is an increasing 
function of Vsd to account for bias-induced decoherence 
effects [SHHZQ], with F -- 14d/log^(V;d/7if) for V^d > 
Tk- 

The expression of the imaginary part of the self-energy 
respects the (particle-hole) symmetry Vsd —t —Vsd, as the 
total dissipation in the system is the sum of the dissipa¬ 
tion coming from electronic and hole processes. We also 
remark that Smll^ becomes substantial for F ~ Wg. In 
this regime, and for eg ~ F, this even function goes to 
zero rapidly for \ Vsd\ > 2F. This behaviour is responsible 
for the dips that occur for Vsd ~ ^ TTr- in Fig. 4, 

as the dissipation steming from the electronic system de¬ 
creases. For I Vsd I > Tk, F is an increasing function of 
Vsd, allowing Q'TOn'^(a;g) to remain substantial, leading 
to a robust tt phase in the reflected signal. For bias volt¬ 
ages close to Ec ~ 0.4meV the model ceases to be valid 
because other energy levels such that (1,1) or (0, 0) can 
be involved. Therefore, in the theory, we include the fact 
that F saturates to a large value close to the ultra-violet 
cutoff Ec- The Kondo resonance has almost collapsed 
due to a very short life-time for Vsd Ec- 

For an 5'[/(4) Kondo effect, one predicts eg ~ Tk and 
for an SU (2) Kondo effect (which is equivalent to a bro¬ 
ken S'C7(4) symmetry onto a lower SU{2) spin symme¬ 
try), one gets eg = 0 [29l[80]. Flere, the theory fits the 
experimental results for eg = Tk (Fig. 4). 

As shown in Figs. 4 and 5, this analysis successfully 
corroborates the experimental results at point A. For 
small bias voltages, the photon field is sensitive to the 
Kondo physics producing Kg Ki and a phase shift close 
to TT in the microwave reflected signal. Interestingly, the 
decoherence effects induced by the bias voltage enlarge 
the value of F for Vsd > Tk, which corresponds to pro¬ 
gressively destroying the Kondo resonance; this is suffi¬ 
cient to maintain the tt phase shift until large values of 
Vsd of the order of the charging energy. Note that the 
light signal is very sensitive to the Kondo effect, as uJq 
and Tk are of the same order of magnitude. Note that 
the low bias features are not present in Fig. 3(c) due to 
the small number of points. Fig. 3(c) only presents the 
general evolution of and Ki at a large energy scale 
{Ec). 

We stress that the datas cannot be explained by a 
resonant-level model, for which F does not depend on 
Vsd- From our impurity model, we indeed know that the 
phase shift from tt to 0 occurs for a bias voltage such 


that eVsd — F. Then the results of Fig. 5(a) would sug¬ 
gest a constant F ~ Ac in the case of a resonant-level 
model. But such a value of F is in contradiction with the 
Coulomb diamond: we would not have such a resolution 
if it were the case. This is also in contradiction with the 
presence of two distinct energy scales in the amplitude 
of the reflected microwave signal. We remark that there 
exists a low bias anomaly concerning the phase response 
in Fig. 5a, which is not present theoretically. From the 
input-output approach we expect that the phase remains 
fixed and equal to tt at low bias. Near this point, the re¬ 
flection coefficient is indeed expected to be real and neg¬ 
ative. The Fermi liquid corrections [45] . which are not 
taken into account in the input-output approach, may 
affect the behavior in this region. 

It is also interesting to investigate the behavior of the 
phase as a function of temperature as shown in Fig. 8. 
One concludes that the phase also allows to probe the 
Kondo temperature and approaches tt for Vsd — t 0 when 
T < Tk- a more precise modelling with both the bias 
voltage and the temperature remains an open question. 
The phase behavior in the regime V ~ Tk is very sensi¬ 
tive to the precise modelisation of the bias and tempera¬ 
ture dependence of F, as also observed in the experiment 
of Ref. [ 3 ] . 


III. COULOMB BLOCKADE REGIME AND 
INPUT POWER DEPENDENCE 


We now discuss amplitude and phase response away 
from charge degeneracy, at point D well inside the 
Coulomb diamond (see Figure 2(a)). Here we have Ki = 
1A32MHZ, Ke = 0.721MHz and = 7.35076GiFz at 
Vsd = 0 and small input power. By increasing Vsd to 
0.8mU, we obtain Ki = 0.564M7Fz, Ke = 1.358MHz and 
Wg = 7.3508GiFz, from a fit result. Figures 9(a) and (b) 
show the phase and amplitude response at a fixed driving 
frequency w = Wg at point D. A phase shift is confirmed 
from bias voltage Vsd = 0 to Vsd — O-SmV. The ampli¬ 
tude sharply goes down for about —15dB from Vsd = 0 
to Vsd — iO.lmU, and two side peaks appear around 
Vsd — ±0.15mU. Then, it continues to go down and two 
dips appear around Vsd — ±0.4toU, which are again con¬ 
sistent with the charging energy Aj, of the quantum dots 
measured by the Coulomb diamond. At small Vsd, this 
confirms that the Coulomb blockade phenomenon pins 
the charge fluctuations on the DQD and as a result the 
photon field is weakly affected by the matter, meaning 
Ke Ki- When the bias voltage increases and compen¬ 
sates the charging energy, there is a revival of charge 
fluctuations on the DQD then resulting in Ke^ Ki. This 
explains the jump of the phase from (almost) zero to 
TT at small bias voltages when increasing the microwave 
power. It is important to note that the cavity is only cou¬ 
pled to the orbital degree of freedom, leading to a phase 
of zero when Vsd —>-0. As a result, it is not sensitive to 
an eventual spin Kondo effect where the spin and orbital 







FIG. 9. (color online). Amplitude (a) and phase (b) response 
as functions of the bias voltage for the point D. Same results 
when increasing the microwave power. 

degrees of freedom are decoupled (resulting from the sys¬ 
tem composed of the left dot and the left lead alone). A 
more refined theoretical analysis would be necessary to 
describe quantitatively the crossover with voltage and to 
include the charging energy effects. 

We then study the power (P) dependence of the res¬ 
onator response as a function of the bias voltage in the 
Coulomb blockade regime at the point D. Fixing the driv¬ 
ing frequency at resonance w = Wq, we show the ampli¬ 
tude response (figures 9(a) and (c)) and the phase re¬ 
sponse (figures 9(b) and (d)) for two different values of 
the power. We see that the phase response is in partic¬ 
ular drastically modified as the power changes. At high 
power the amplitude and phase response are close to the 
ones obtained at point A. This observation can be ex¬ 
plained by the fact that the photonic interaction acts as 
an artificial chemical potential and may compensate the 
charging at high drive powers. Then the high drive power 
may bring the system from the Coulomb blockade regime 
to the resonant regime. 

At a mean field level, the equivalent charging energy 
reads Ec—X{a+a^). In the steady state, the input-output 
theory allows to relate the drive power to the number of 
photons in the cavity and ultimately to (a -I- a^) if we 
suppose that the steady state of the cavity is a coher¬ 
ent state. Following Ref. we have P = hujj{a^a), 
where t] is the frequency-independent decay rate of the 
photons in the leads. For a coheren t state in the cav¬ 
ity this allows to write (a'^ -|- a) ~ ^P/{hujr]). We can 
then estimate the critical power P^. needed to compen¬ 
sate Ec- We have Ec = OAmeV, w = 7GHz, and 
we take A ~ 0.15GHz and 77 ~ O.lGHz. This finally 
gives Pc — 10~^^V = —70dBm, which is consistent 
with the experimental observations. We note that the 


high-power responses at small biases differ slightly from 
those at point A. This may be linked to additional non¬ 
equilibrium effects arising at high power. 

We provide additional theoretical derivations in Ap¬ 
pendix B, where we also check that bringing the system 
from a Coulomb blockade regime to a resonant regime 
is not possible for the points B and C. In this case, the 
system always remains in a blockade regime as a result 
of the strong applied orbital magnetic field. 


IV. CONCLUSION 

In conclusion, we have measured the response of a su¬ 
perconducting resonator capacitively coupled to a biased 
graphene double quantum dot (DQD). This setup allows 
to probe dynamically the electronic degrees of freedom 
with light. We report a full spectroscopic analysis of the 
resonator’s response as a function of the DQD bias and 
gate voltages. We have furthermore studied the power 
dependence in the phase and amplitude responses at var¬ 
ious points. We explain the experimental findings with 
a Kondo impurity model entangling the electronic spin 
and orbital degrees of freedom m taking the graphene 
band structure into account. We assumed that the chem¬ 
ical potential does not lie at the charge neutrality point 
in order to have a finite density of states in the reser¬ 
voir leads, making the Kondo energy scale Tk finite. We 
emphasize that the signatures of such a Kondo effect is 
visible, both in the phase tt at low temperatures and in 
the dips in the amplitude of the reflected signal occurring 
roughly at Vsd ~ ETk- This is also consistent with DC 
transport datas across the DQD, exhibiting a small zero- 
bias anomaly at Vsd which reflects the fact that 

the Kondo resonance is peaked above the Fermi energy 
in this system. Note also that the phase of tt in the re¬ 
flected signal is distinct from the phase of 7r/2 predicted 
by the theory in the transmission of light [55]. Light 
spectroscopy can then probe the dynamical properties of 
correlated electron phenomena when the light frequency 
is (almost) synchronized with the typical frequency scale 
associated with the correlated phenomenon. Our results 
offer a platform for tunable hybrid quantum electrody¬ 
namics and non-equilibrium impurity physics. 
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Appendix A: General Arguments on the Light-Matter Coupling 

The general Hamiltonian of the double dot setup (at any point of the charge stability diagram) takes the form [50] 
Hi = Hl + Hu + Hd + Ht, where 


Hi=L,R = + Vl)c\^iCkal 

k<7 

Hd= '^{^pid^piadpia) +lLnL{nL-^)+lRnR{nR-l)+^MnLnR 

l—L,R picr 

Ht= ^ '^'^{Hnicl^ldp^a + h.C.) + ^ {'^PlPr dLad pj^a + h.C.). (-^ 1 ) 

l—L,R ka PI PLPRiy 

The first line describes the electrons in the Left (L) and Right (R) leads, and Ckai (Ck^i) is the annihilation (creation) 
operator of an electron of spin cr in the mode k and in the lead 1. The second line is the Hamiltonian of the two dots: 
dp^. «.) is the annihilation (creation) operator of an electron of spin a on the dot I = L,R and at the energy level 
pi- Here the y’s are functions of the capacitances of the tunnel junctions and ni denotes the number of electrons in 
the dot 1. The last line accounts for lead-dot and inter-dots tunneling events. We now take into account the coupling 
of the left lead to the superconducting resonator, so that the total Hamiltonian becomes H = Hi + H 2 with 


H 2 = ujqo^ a + \{a + a^) I 

\ k<7 


^kaL^kaL 


We apply one first unitary transformation H' = W HU, with 


U = exp 


^kcrL^^kaL 
\ k(7 J 


(A2) 


(A3) 


where (j) = ° ^ . We can compute the transformed Hamiltonian H' thanks to the Baker-Campbell-Haussdorff 

formula e^Ye~^ =Y + [X, Y] + [X, [X, Hj] -|- [X, [X, [X, H]]] -|- ... We then get 


U^CkalU = Ckaie 


i(f) 


{uJoa)a) U = -A(a + a'^) CkaL<^ka-L^ 

A(a + a^) E ^kaL^kuL 




k<7 




(A4) 

(A5) 

(A6) 


Assuming small flucutations concerning the total number of conduction electrons, the right hand side of Eq. (6) is 
almost a constant, so that the resulting Hamiltonian reads H' = Hl + Hr + Hr + H'rp + wga^a, where we have 


^t = EE(^ kpL^kak'^PL'y^ + y]] + ^-C.) + {^PlPr d^PR^d Pfi<T T h.C.). (A7) 

k<7 PL ka PR PLPRa 


This unitary transformation has suppressed the explicit coupling between the left lead and the resonator. The 
tunneling terms from the Left lead to the Left dot have moreover acquired a phase which depends on the state of the 
resonator. 

We apply then a second unitary transformation H = V^H'V, with 

( E d-lLadpLa 

\pL<y 


V = exp 


(AS) 
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After similar calculations, the resulting Hamiltonian reads H = Hl + Hr + Hu + Ht + where we have 

Hu = HD-Ka + a^) 

\pL<y / 

y^jtkpicl^idpicr + h.C.) + ^ 

i^PLPR^PLCr^PR<^^ )■ (A9) 

l=L,R ka Pi PlPr<^ 

In this final form, the mode of the resonator is coupled to the energy levels of the Left dot. Additionally, the 
tunneling term between the dots has acquired a phase which depends on the p operator of the mode. 

This analysis implies that the form of the reflection coefficient in Eq. (5) of the main text is valid for all the 
points in the charge stability diagram. 


Appendix B: High power regime 


1. General arguments 

For large drive, we decompose the photon annihilation operator as a = (a) + d. The hrst part (a) is determined 
by the strength of the drive field and the damping rate, and we suppose that steady state of the cavity is a coherent 
state. We see that the constant semi-classical part of the drive (a -I- a^) acts as a chemical potential on the Left dot. 
Increasing the drive power will then move the system on the honeycomb phase diagram of the double dot. Following 
Ref. [74] the chemical potentials of the two dots pL and pu, read as a function of the number of electrons on the dots 
{Nl and Nr), 


— X{a + Ecl + NrEc^ — {CghYLpEc^ + CgRVRpEc„^) (Bl) 

Nr) = ^Nr - Ecu + NrEc,^ - {CgLVLpEc„, + CgpVRpEcji ), (B2) 

where Ecj is the charging energy of the individual dot j, is the electrostatic coupling energy, Cj (Cgj) is the 
capacitance coupling the dot j to the neighboring lead (gate), and Cm is the capacitance which couples the dots 
together. The Left and Right gate voltages are denoted by Vpp and Vrr. The effect of the driving can be absorbed 
into the definition of new gates voltages V[p and Vrr, 


Pl{Nli Nr) = 




V[p = Vlp 
^RP = ^RP 


X{a + a'^)|e| 

C,L (Ec, - gj) 

X{a + a^)|e| 

CgR [Ec^ - ■ 


(B3) 

(B4) 


Driving the system at high power allows to move the state of the system on the honeycomb phase diagram along a 
line characterized by the equations (|B3| and (|B4|). This corresponds to 


SVrp _ CgL Eg^ 
SVrp CgR Ecf^ 


(B5) 


Such a line has a slope which is the opposite of the line of degeneracy between the states {N, M) and {N, M + 1) 
(the lines where lies the point B). 


2. Theory at Point A 

We can see from Fig. 10 that driving the cavity may bring the system at resonance (close to the point A). 

Now we give an estimate for the input power that allows to make this transition from the Coulomb Blockade regime 
to the resonant regime. In the steady state, the input-output theory allows to relate the drive power to the number of 
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Vrp (V) 


FIG. 10. (color online). Effect of the driving on the state of the system. The full red line shows the shift in the phase diagram 
induced by the drive at point D. The dotted red lines accounts for the shifts caused by the drive from the points B and C. 


Point B 


1 -0.5 0 0.5 1 

^'sd (mV) 

Point B 

/ \ — -TOdBm 


-1 -0.5 0 0,5 1 

(mV) 





FIG. 11. (color online). Amplitude (a) and (e) and phase (b) and (f) responses as functions of the bias voltage for the point B 
and G. Amplitude (c) and (g) and phase (d) and (h) responses when increasing the microwave power for the point B and G. 


photons in the cavity and ultimately to (a + a'^). Following Eq. (E45) of Ref. [75], we have P = {huji]/i) {a'^a), where 
r] is the decay rate of the photons in the leads (assumed to be independent over the range of frequencies relevant to 
the cavity - Markov approximation). If we suppose that the steady state of the cavity is a coherent state, we have 


{a^ + a) ~ 



(B6) 


We can then estimate the critical power Pc needed to compensate Ec- We have Ec = 0.4meE, uj = IGElz^ and we 
take A ~ O.lSGi/z and ry ~ 0.1G77z. This finally gives 


Pc ~ —70dBm. 


(B7) 


This result is consistent with the experimental observations. 


3. Experimental Checks at points B and C 


Bringing the system from a Coulomb blockade regime to a resonant regime is not possible for the points B and 
C, as illustrated by the dotted red lines in Eig. 10. The system always stays in a blockade regime. This behavior is 
verified experimentally, as shown in Eig. 11. Driving the system at high power does not affect (much) the response 
of the reflection coefficient. 
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